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Abstract 

We argue on a basis of a simple few mode model of SU(2) Yang-Mills theory 
that the color off-diagonal coupling of the soft plasmon to hard thermal exci- 
tations of the gauge field drives the collective plasma oscillations into chaotic 
motion despite the presence of the plasmon mass. 
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The chaoticity of classical nonabelian gauge fields was originally discovered for spatially 
homogeneous gauge potentials [jXJ] and later extended to special, radially symmetric field 
configurations 0,0]. More recently, the chaotic dynamics of the full classical Yang-Mills 
equations was investigated in numerical studies of Hamiltonian SU(iV c ) lattice gauge theory 
with N c = 2, 3 color degrees of freedom The complete Lyapunov spectrum of the SU(2) 
gauge field was calculated || , and it was found that the field configurations ergodically cover 
the available phase space on the time scale of the inverse of the leading Lyapunov exponent, 
r pa 1/Ao oc l/g 2 E, where E is the energy per lattice degree of freedom. The leading 
Lyapunov exponent coincides with twice the plasmon damping rate obtained in resummed 
perturbation theory at finite temperature T 0. 

The coupling to the heat bath gives rise to a mean inertia of the plasmon field. For 
long wavelength {k — > 0) and static fields (ou = 0) it coincides with the Debye mass square 
mp = ^g 2 T 2 ; for on-shell plasmons in the limit k — > one has mp = ^-g 2 T 2 . Such a mass 
restricts the region of chaotic motion to large amplitude plasmon fields. In fact a control 
parameter, 

with E being the total classical energy of the plasmon field and V the plasma volume, has 
been found to control chaos in simple, few-mode models of the Yang-Mills field theory ||. 
For small R < 0.2 chaos develops. If the mean plasmon mass mp were the sole effect of hard 
thermal gluons, then small amplitude plasmon fields would not grow chaotic |J. 

We now argue that even small amplitude collective plasma oscillations can be chaotic 
because of the color-nondiagonal coupling to the fluctuating components of the plasmon self- 
energy. These color non-diagonal matrix elements can render some eigenvalues of the mass 
square matrix Mp negative. This is a source of instability for small amplitude soft fields. 
Once the energy stored in the soft components has grown, the classical chaotic behavior sets 
in. 

We are not able yet to prove the above scenario generally. In this paper we study as an 



example the coupled dynamics of some selected Fourier components of the vector potential 
to show how the mechanism discussed above can be realized. The central point here is the 
existence of a hierarchy of relevant timescales in the dynamics of oscillations of the gauge 
field. In fact, three such timescales can be identified.^ The hard thermal modes oscillate 
over a period of t\ ~ 1/T. Due to the nonabelian coupling they influence the soft mode 
dynamics on the timescale £2 ~ I/9T. This also sets the timescale for the evolution of the 
soft plasma oscillations. Finally, the thermal ensemble average over particular hard mode 
configurations is characterized by the timescale of the gluon damping rate £3 ~ l/g 2 T. The 
fact that £2 < £3 requires that the ensemble average has to be taken over solutions of the 
plasmon equation of motion rather than over the equation itself. 

Let us now consider a simplified model retaining the essential features. Here we build on 
the infrared limit of Yang-Mills field theory ( "Yang-Mills mechanics" |IJ). Each Fourier mode 
of the SU(2) gauge field representing hard thermal gluons carries six degrees of freedom (two 
polarizations, three internal degrees of freedom). Soft modes, neglecting spatial derivatives, 
can be described by the three eigenvalues of the 0(3) x 0(3) — symmetric vector potential 
Af. Without loss of generality they can be chosen as [[IT]]: 

x = A{, y = A 2 2) z = A\. (2) 

The most general ansatz describing one hard and one soft momentum contains therefore 
nine coupled modes 

A\ = w\, A\ = 0, 

A 2 2 = y + wl Al = 0, 
Al = wl Al = z. (3) 

Here 



A\ = x + w{, 
Al = wl 
Al = wl 



1 We here assume, as is customary in the analysis of thermal gauge theories, that g is sufficiently 
small that the hierarchy t\ <C £2 ^ £3 holds. 
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wf = y/2 (v,i cos(k ■ r) +u?sin(fc-f)) (4) 
denote the hard thermal modes with time-dependent amplitudes u\ and i>". We choose 

— * 

the thermal wave number fc = Te 3 to point in the third direction. Then all longitudinal 
components w% = vanish and only one derivative (in the third direction) differs from zero 
in this ansatz: 

dawf = TV2 (-< sin(£ • f) + vf cos(£ ■ f)) . (5) 

Choosing the temporal gauge, A% = 0, the electric field components are simply the time 
derivatives of the vector potential 

E\ = x + w\, El = w\, E\ = 0, 

El = wl E% = y + wl E 2 3 =0, 

E\ = w\, El=wl El = z. (6) 

The nonabelian magnetic field has the following components 

B\ = -d 3 wl + gz(y + wl), 

B\ = d 3 w[ - gzw\, 

B\ = g(w\wl - wlw\) - gywf, 

B\ = -d 3 wl - gzw\, 

B\ = d 3 wl + gz(x + w{), 

B l = 9{w\w\ - wlw\) - gxw\, 

B\ = -d 3 wl 

Bl = d 3 wl 

Bl = g(w\wl - w l 2 w\) + g(w{y + w 2 2 x) + gxy. (7) 

In order to obtain the effective Hamiltonian for the soft modes, which are taken to be 
constant throughout space, we perform a spatial average over the relevant terms in the 



Hamiltonian. For the quadratic and quartic expressions of the hard amplitudes wf we use 
the following relations 



wfdsw^T^-vtuf), 
dswfd^ = T 2 (ulu) + vtvfj , 

wfwfafw* = - (u»5«J + vyy^ + «X u ? u j + u »j v j) ■ (8) 



The spatially averaged Hamiltonian, 



^ 1- 



= + -Bf Bf ^Ht + Hs + Hj (9) 

contains three parts: ify, describing the hard thermal modes interacting among themselves; 
H s , describing the soft modes in isolation; and Hj denoting the interactions between hard 
and soft modes. According to the usual picture of a gluon plasma we treat the hard modes 
perturbatively, neglecting their mutual interactions. This yields: 

H T — > H = \ £ £ [(« + «X) + T 2 (« + «)] , (10) 

Z i=l a =l 

in the leading order approximation. 

The part of the Hamiltonian which contains soft modes must be treated non- 
perturbatively, so we keep the quartic anharmonic terms here. We get 



H S = \ (f + f + z 2 ) + 9 - (x 2 y 2 + y 2 z 2 + z 2 x 2 ) (11) 

for the soft mode part and 

H t = ^ (Mix 2 + M 2 y 2 + (Ml + M\)xy + M 3 V) + Q 3 z (12) 

for the interaction between soft and hard thermal modes. Here 

Ml = g 2 {u\u\ + vlvl + u\u\ + v 2 2 vj) , 
M l = 9 2 ( u lul + vfvf + u\u\ + vlvl) , 
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u\u 2 — v it; 2 



2"i "2^1 y > 



M 2 = Ml = g 2 (2u\u\ + 2v\i 
Ml = g 2 (u]u\ + v\v\ + u\u\ + v\v\ + u\u\ + v\v\ + -u^u 2 , + f 2 ^) > 
g 3 = 2#T f it Jf 2 - u\v\ + ^ - ufo 1 ) (13) 



are the components of the effective mass square matrix for the soft modes and Q 3 is the 
fluctuating color charge of the hard modes, which acts as a source for the soft modes. 

The Hamiltonian Ho (|T0| ) describes independent oscillations with a frequency of uo = 
\k\ = T. For the dynamical evolution of the soft modes x, y, z this is a very fast motion; we 
can therefore average over many time periods in order to arrive at an effective soft mode 
Hamiltonian Hi for the interaction. The ensemble averaging over the relative phases and 
slowly fluctuating amplitudes of the hard thermal components, however, cannot be done at 
this point, because these quantities change on the timescale t 3 ~ l/g 2 T. Consider the hard 



thermal eigenmodes of the Hamiltonian (10) 



<(f) = E?cos(Tf + tf), 

vt(t) = V? C os(Tt + r i ). (14) 

The initial amplitudes Uf, V°" and phases ip®, ipf (i — 1, 2 a — 1, 2, 3) determine the quantities 
Mf and Q a which in turn influence the evolution of the soft modes x, y, z. In order to obtain 
the correct scaling Mf ~ g 2 T 2 of the elements of the mass matrix flT3"| ) we take the amplitudes 
of the fast oscillations to be 0(T). As seen from (TTO) this is consistent with an energy density 



H ~ T 4 in the hard modes, as expected for the thermal bath. 

If, for the sake of simplicity, we assume all amplitudes equal to each other, so Uf = V°" = 
U ~ T, we obtain M 1 = Mf = m 2 , Mf = 2m 2 with 

m 2 = 2g 2 U 2 ~ g 2 T 2 . (15) 

These effective mass terms in the averaged Hamiltonian Hj are the analogue of the plasmon 
mass term in the full thermal gauge theory. For small amplitudes of the soft modes x, y, z 
they dominate over the anharmonic terms in ( |TTD and apparently cause the soft modes to 
oscillate harmonically with a frequency u>s ~ gT. 
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However, the presence of non-diagonal elements Mf in the mass square matrix can desta- 
bilize these oscillations. The time averaged off-diagonal term Ml, as well as the source term 
Q 3 , depend on the phases of the hard thermal oscillations (|T4|): 

m 2 / 1 1 \ 

Ml = — (cos(<^ - <pl) + cos^ 1 - ^2) - 2 cos (^ - tf) ~ 5 cos (^2 - A 2 )J , (16) 

and 

Q 3 = ~y~ ( cos ifl ~ - cos(<^2 - V2) + cos(^ - - cos(v?i - ^i 1 )) • (17) 

Because of their dependence on the relative phases of the fast modes Ml and Q 3 would 
vanish if a complete ensemble average over the fast modes were performed. However, as 
pointed out at the beginning of our discussion, the phases and amplitudes of these modes 
change more slowly (on timescale £3) than the soft modes oscillate (on timescale £2)- Hence, 
both Ml and Q 3 must be considered as secular quantities in the averaged (over timescale 
ti) Hamiltonian H$ + Hj governing the evolution of the soft modes. 

Far from the "average" case, Ml = Q 3 = 0, there is the extreme possibility of M\ = 
3m 2 /2 and Q 3 = m 2 T/2g ~ gT 3 obtained for <p\ = <p% = ipf = tpl±ir, ip\ = ^)\ = ^\ = <^±vr. 
As we will see below, in this case the mass square matrix of the coupled x and y modes 
has negative eigenvalues giving rise to an exponentially growing solution of the linearized 
equations of motion for small amplitudes x and y. 

We now consider small amplitude oscillations and neglect terms nonlinear in x or y. The 



z-coordinate then decouples from the x — y motion, as seen from (0). It is most convenient 
to cast the resulting equation into matrix form 



% + Ml MA (x\ (fr 



(18) 



\Ml ^ + Ml)\y) V /' 
where we remind the reader that the "bar" denotes the time averaged quantities for fixed 
phase <fi,ipi and amplitudes U^Vf 1 . Since the average over the ensemble of these initial 
conditions must be performed after the full solution is obtained, we proceed with solving 
(|18D by finding the eigenvalues and eigenvectors of the matrix operator on the left hand side. 



Seeking the solution with x, y oc e luJt , we find the eigenvalue equation 

u A -u 2 (Ml + Ml) + (MlM$-M?M£) = 0. (19) 

Using the previous result M\ = Mf = m 2 and = M\ = r ym 2 with |7| < 3/2, we find 

^ = m 2 (l±| 7 |). (20) 

Clearly, for \^\ > 1 one of the eigenvalues becomes unstable with an imaginary frequency. 

With the help of the eigenvalues uj± and eigenvectors (x±, y±) the general solution of the 
time-averaged equation (|T3| ) can be written in the form 



where 



. , s ^ sinc^t / c 2 adA 

A(t)= £— U <?)■ (22) 

where (q, dj) are the eigenvectors corresponding to the eigenvalues w?. It is not difficult to 
show that the product Qcfj vanishes when it is averaged over the phases tpl and As a 
result, the full thermal ensemble average gives rise to a color diagonal effective soft mode 
propagator (A(£)), which contains an exponentially growing part originating from those 
regions of phase space where lo 2 _ < 0, as is seen from (^). 

Summarizing, we demonstrated in a simple model that the coupling of soft plasmon 
oscillations to hard thermal gluons can drive the soft field amplitude into an exponential 
growth and with that into a chaotic dynamical behavior. There are, of course, many possible 
few mode approximations to the Yang-Mills field theory. The example we discussed above 
was generic in the sense that it showed both dynamical mass generation and the development 
of chaotic instability due to off-diagonal parts of the self-energy, as well as source term 
coupling to the soft modes. One may expect that by this mechanism the dynamics of long 
wavelength modes in a hot gluon plasma remains chaotic despite the presence of a dynamical 
plasmon mass. 
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